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Within the t-J model Hamiltonian we present a RVB mean field theory directly in terms of dopant 
particles. We apply this theory to Naa;Co02 ■ y^20 and show that the resulting phase diagram Tc 
versus doping is in qualitative agreement with the experimental results. 



I. INTRODUCTION 

The remarkable discovery of superconductivity in 
Na^Co02 • 2/H2O for x = 0.35 and y = 1.30 by Takada 
et al} attracted a lot of attention. The experimental 
findings indicate several striking similarities between the 
cobaltates and the cuprates. Naa;Co02 • can be 

viewed as a 2D Mott insulator. The Co atoms form a tri- 
angular lattice but Co**^ is in a s = ^ low spin state. The 
transition temperature {Tc) is seen to decrease for both 
underdoped and overdoped materials^ although for the 
cobaltates the maximum Tc is much lower {Tc « 5if ) and 
the optimal doping is twice as large as in the cuprates. 
Finally, as one varies temperature and electron concen- 
tration, apart from superconductivity, there are observed 
unusual electronic properties'^ and clear hints that strong 
electronic correlation is at work in both cases. 

All this turns attractive the application of the res- 
onating valence bond (RVB) ideas to this new com- 
pound. Baskarani was the first to present qualitative 
arguments in favor of the RVB approach for the cobal- 
tates. Soon after that Kumar and Shastry- as well as 
Lee and coworkers?'^ presented their first estimates for 
the mean field (MP) phase diagram, Tc versus doping, in 
a RVB framework. In the cobaltates the C0O2 layers are 
arranged in a triangular lattice which naturally exhibits 
considerable magnetic frustration. This not only brings 
modifications to the symmetry of the resulting supercon- 
ducting state, as pointed out by otherS )'^i^i^i''i^ but makes 
the the application of a MF RVB, to reproduce the ex- 
perimental findings related to the superconductivity of 
Naa;Co02 •yH20, even more challenging. The reasons for 
that are as follows. The maximum Tc for the C0O2 lay- 
ers occurs for doping values nearly twice as large as in 
the cuprates. This might be indicative that the phase 
fluctuations of the superconducting (SC) order param- 
eter could be much too strong for the stability of the 
MF RVB state. Moreover, in the standard approach low 
doping is always favored and makes even harder a more 



quantitative agreement with experiment in the case of 
the cobaltates. 

Within the original Baskaran-Zou- Anderson (BZA) 
MF approximation^ a non-zero value of the RVB MF 
order parameter (OP) does not by itself imply supercon- 
ductivity. The true SC OP in their approach is essen- 
tially taken as a product of a spinon (a spin-1/2 neutral 
fermion) pairing OP and a bose condensation factor for 
the holons (spin-0 charged bosons) fiS following the more 
conventional slave-boson approximation. The phase of 
the OP accounts for the fluctuations which drives Tc to 
zero at zero doping. The bose condensation temperature 
for the holons is estimated separately and the region in 
which both the spin-pairing OP and the mentioned holon 
bose factor are non-zero determines the resulting RVB 
superconducting phase. ^ All this seems to indicate that 
the standard MF RVB basic ingredients - the lattice spin 
singlet pairs - which might indeed be appropriate to de- 
scribe the physics at very low dopings is not the best 
starting point to address the superconducting regime at 
higher doping values. 

To overcome those limitations we present a RVB MF 
scheme which takes direct account of the dopant par- 
ticles themselves and treats the non-double occupancy 
(NDO) constraint beyond the conventional slave-boson 
mean field approximation. As will be demonstrated by 
this and later works this RVB representation is most suit- 
able to deal with cases in which strongly correlated elec- 
tronic superconductivity is manifest in both low and high 
doping regimes. 

In this first work we apply our method to the C0O2 
superconductors. Our starting point is the t-J model 
on a triangular lattice. In doing that we follow the argu- 
ments which consider the 3d levels of the Co''^ ions being 
crystal field split in the C0O2 layers producing singly oc- 
cupied non-degenerate spin-1/2 orbitals. Those or- 
bitals are directly associated with the singlet states in 
our t-J model representation based on the Hubbard X 
operators. The use of that representation will allow us to 
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go beyond the conventional treatment of the NDO con- 
straint. The X operators are later given in a convenient 
coherent-state path-integral representation. The result- 
ing variables in this representation are naturally split into 
bosonic and fcrmionic degrees of freedom. The bosonic 
modes correspond to SU(2) spin excitations while the 
fermion variables are spinless and describe U(l) charge 
excitations instead. Combining those spinon amplitudes 
and spinless fermion parameters together we then con- 
struct appropriate fermionic fields which carry both spin 
and charge degrees of freedom and can be directly re- 
lated to the dopant carriers in the t-J model. We are 
able in this way to take a direct account of the doping 
dependence of the critical superconducting temperature 
preserving all the symmetry properties of the t-J Hamil- 
tonian. 

We reformulate the RVB theory of the SC phase en- 
tirely in terms of those quasiparticle states and use this 
scheme initially to describe the superconducting proper- 
ties of the cobaltates. We find qualitatively good agree- 
ment with experiment and we are able to reproduce the 
observed dome structure of the versus doping phase 
diagram for those materials in a RVB MF framework. 
The application of our RVB method to the cuprates will 
be presented in a subsequent work. 



II. t-J HAMILTONIAN AND THE NDO 
CONSTRAINT 

We start by expressing the t-J Hamiltonianii 
Ht-j = -tY^cl^Cj^ + h.c. + f^jt^j - ^n^rijj , 

ija ij 

(1) 

with the NDO constraint, ^„ni„ < 1, in terms of the 
Hubbard operators/^ 

= (1 - "j-<t) , rii-^n^a- = 0. 

Here Cia is the electron annihilation operator at site i 
with the spin projection a =11, Uia- = cf^Ci^, and the 
Q's are the corresponding electron spin operators. In 
terms of these operators the local NDO constraint holds 
rigorously and the t-J model becomes 

ija ij 

where the electron spin operator now reads Qi = 
^T,c7c7' ^i°^<y<y'^i''' ^ with the r*'s being Pauh matri- 
ces. 

Fermionic operators X^'^ project the electron creation 
operators onto a space spanned by the basis {|0)i, |cr)i} 
and take the form = \(t).,{Q[,. Together with the 

bosonic generators, X^'^ = \(T)i{a'\i the full set of op- 
erators X""^, a,b = 0,t,i forms, on every lattice site, a 



basis of the fundamental representation of the semisim- 
ple doubly graded Lie algebra su(2|l) given by the 
(anti)commutation relations 

{Xf,Xf}± = {XfS'" ± (5"'*)(5*^ 

where the (-I-) sign should be used only when both oper- 
ators are fermionic. 

Since su(2|l) can be viewed as a supergeneralization 
of the conventional spin su(2) algebra, the t-J Hamil- 
tonian appears as a superextension of the Heisenberg 
magnetic Hamiltonian, with a hole being a superpart- 
ner of a su(2) magnetic excitation^ This superalge- 
bra can also be thought of as a natural generalization 
of the standard fermionic algebra spanned by genera- 
tors c+,Cct, and unity /, to the case where the fermionic 
operators are subject to the NDO constraint. The in- 
corporation of this constraint manifests itself in more 
complicated commutation relations between X operators 
in comparison with those produced by the conventional 
fermionic operators. Note that the Gutzwiller projection 
Pg = Yli (1 ^ nia-rii^cy) that excludes the doubly occu- 
pied state I ti) is equivalent to the Hubbard operator 
representation, since PcclyCjaPG — X^^X^" . 

Note also that the occupation constraint is different 
for the hole and electron doping. To treat them in a 
unique way we perform, for electron dopings, a canonical 
particle-hole transformation Cia cf-a that restores the 
non-double occupancy constraint but reverses the sign 
of t. Using then the Hubbard operator representation 
in terms of the transformed c-operators we again arrive 
at Eq.Q with, however, t — > ~t. Although the C0O2 
layer is an electron doped Mott insulator we shall for 
convenience formally deal with the more familiar case of 
hole doping making the necessary changes only at the 
end of our work. 

Since the X operators are generators of the su(2|l) su- 
peralgebra we are lead naturally to employ the su(2|l) 
coherent-state path-integral representation of the t-J 
partition function. There are a few rationales to do that. 
First, this provides a mathematical setting well adjusted 
to address the t-J model with the crucial NDO constraint 
naturally built in the formalism from the very beginning. 
Second, within the su(2|l) path-integral representation 
the associated effective t-J action lives on a natural clas- 
sical phase space of the t-J model - the SU(2|1) homoge- 
neous compact manifold, CP-'^I^ (see below). The group 
SU(2|1) acts on the CP^'^ manifold as a group of canon- 
ical transformations in a way that the transformation 
properties of the basic fields - the local coordinates on 
CP^I^ ~ can be easily found. Third, these coordinates 
are naturally split into bosonic and fermionic degrees of 
freedom. In the context of the t-J model the bosonic 
fields correspond to the SU(2) spin excitations whereas 
the fermionic ones are spinless and may be used to de- 
scribe the U(l) charge excitations. This provides a natu- 
ral setting to implement the spin-charge separation inher- 
ent in the spin liquid phase at least in ID. Finally, the 
transformation properties of the CP^'^ coordinates un- 
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der global SU(2)x U(l) rotations - the exact symmetry 
of the t-J Hamiltonian - imply that their certain combi- 
nations transform in the linear spinor representations of 
SU(2) and may therefore be used to describe fermionic 
quasiparticle excitations that carry both the charge and 
spin quantum numbers. Wc show that such quasiparti- 
cles arise as the dopant particles in the t-J model relevant 
for describing the SC phase. In particular, we formulate 
the RVB theory of the SC phase of the t-J model directly 
in terms of the dopant particles and apply it to describe 
the SC in the cobaltates. 

Some earlier attempts have also been made to apply su- 
persymmetry to study many-fermion interacting Hamil- 
tonians as well as the t-J and related models. A lin- 
earization scheme for the general Hamiltonian of an in- 
teracting fermion system has been proposed in Ref. 0. 
A hierarchy of spectrum-generating algebras and super- 
algebras including su(2|l) results from such a new mean- 
field treatment. A supersymmetric representation of 
the Hubbard operator which unifies the slave-boson and 
slave-fermion representation into a single U(l|l) gauge 
theory has been developed in Ref. Hsl Such a represen- 
tation makes unnecessary the choice between a bosonic 
and fermionic spin and is most suitable to describe the co- 
existence of strong magnetic correlations within a param- 
agnetic phase. Besides, it has been demonstrated that 
thus defined supersymmeric Hubbard operators prove to 
be very efficient in treating the physics of the infinite U 
Hubbard model. 



where 



d^lsu{2\l) 



dzdz dS^dS^ 



2m l + 



stands for the SU(2|1) invariant measure on the coherent- 
state manifold, CPi|i=SU(2|l)/U(l|l), and / is the iden- 
tity operator in the projected Hilbert space. Explicitly, 
we have 



d^i\zC){z£,\ = 



dzdzd^dS^ 



27ri(l-h|z|2 + eO 
dzdzdS^dS^ 1 



x(i T)(Ti + kni)ai+eeio)(oi) 

= IT)(TI + U)ai + |0)(0|^/. 



In the basis 1^,^) — Hj NjiO)' ^^'^ partition fmic- 
tion takes the form of the su(2|l) CS phase-space path 
integral, 

Zt-j = tr exp(-/3i7t_j) = / D/^sc/(2|i) (2, ^) e^'-' , 

(4) 

where 



D^isu^2\l){z,0 = Jl 



dzj{t)dzj{t) d£,j{t)d£,j{t) 



3,t 



2m 1 + \z,V+Uj 



III. SU(2|1) COHERENT STATES AND PATH 
INTEGRAL 



The normalized su(2|l) coherent state (CS) associated 
with the 3D fundamental representation takes the form 

1^,0 = (l + z-z-HC6-'/'exp(zXiT I (3) 

where z is a complex number, and ^ is a complex Grass- 
mann parameter. The set {z,^) can be thought of as 
local coordinates of a given point on CP^'^. This super- 
manifold appears as a A'^ = 1 superextension of a com- 
plex projective plane, or ordinary two-sphere, CP""^ = 5'^, 
to accommodate one extra complex Grassmann parame- 
terJX At ^ = 0, the su(2|l) CS reduces to the ordinary 
su(2) CS, = 0) = \z) parametrized by a complex 
coordinate z G CP^. Note that the classical phase space 
of the Hubbard operators, CP"'^'^, appears as a A = 1 
superextension of the CS manifold for the su(2) spins. 

As is well known the key point in constructing the 
coherent-state path-integral representation of a partition 
function is the resolution of unity or, equivalently, the 
completeness relation for the coherent states. In terms 
of the normalized set of states © it takes the form 



stands for the SU(2|1) invariant measure with the 
boundary conditions, 2:^(0) = Zj{(3), £,j{0) = —^jiP). 
The t-J effective action on CP^'^ now reads St-j = 
— {z,^\d/dt + Ht^j\z,S,)dt, which gives 



1 ZjZj 



■E 



1 + + 



Hf_jdt. 



(5) 



The first part of the action Q is a purely kinematical 
term that reflects the geometry of the underlying phase 
space while the classical image of the Hamiltonian (0) be- 
comes an average value of Ht-j over the su(2|l) coherent 
states. 



Hf_j = {z,^\Ht-j\z,0 



^i^jil -\- ZjZi) + h.c. 



(6) 



The fact that the electron system with the NDO con- 
straint lives on the compact manifold, supersphere CP^'^ 
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can be explained as follows. Let us for a moment sup- 
pose that the so-called slave-fermion representation for 
the electron operators is used, i.e., 



= f^aL 



(7) 



where fi is a on-site spinless fermionic operator, whereas 
Qia is the spinful boson. The NDO constraint now reads 
l^a'^ta'^io- + ftfi — 1- Within the slave-fermion path 
integral representation 



Zt-J = / Dflfia 



J) 



(8) 



with the integration measure Dfifiat = 
Da.i'jDai'jDai^DaiiDfiDfi , this constraint trans- 
forms into 



fifi 



(9) 



remaining integration we have (the site dependence for 
the moment being suppressed), 

Da^D-a^DfDf - sdet]] "i'/' || DzDzd^D^ 

The Jacobian of the change of the supercoordinates ap- 
pears as a superdeterminant of the transformation ma- 
trix 



.18 



sdetl 



Here 



d{ai,ai,f, /) | 



( dai da I 
Ski 
dz dz 



sdet ^ ^ 
det{A - BD-^C)detD'^ . 



B = 



da I da I 
da I da I 



with Cio- and fi standing now for complex numbers and 
complex Grassmann parameters, respectively. Equation 
is exactly that for the supersphere CP^'^ embedded 
into a flat superspase. Any mean-field treatment of © 
should respect this constraint, which, however, poses a 
severe technical problem. If one however resolves this 
equation explicitly by making the identifications 



f^ = 



1 + ZiZi+ ^j^i 
1 + ZiZi + CiC. 



=7 — 



1 + ZiZi + 



(10) 



one can further treat the variables Zi,S.i as if they were 
indeed free of any constraints. 

Note that the electron operator {Tj) is invariant under 
a local gauge transformation. 



O^ifj^ : fi ^ fi^ ; 



or equivalently, under the change — > (pi + 9i. This 
gauge symmetry is a consequence of the redundancy of 
parameterizing the electron operator in terms of the aux- 
iliary boson/fermion fields. In contrast, the sm(2|1) pro- 
jected coordinates 



Zj=ajx/aj-r, ^i^fi/a, 



are seen to be manifestly gauge invariant. The domain 
of the flat measure in ((SJ that involves the spin up 
bosonic fields can be rewritten at every lattice site as 
Dai^Dtti^ — D\ai^\'^D(f)i. The |ai|p field can easily be 
integrated out from eq.® because of the constraint 
Since the t-J action is U(l) gauge invariant and hence in- 
dependent of (pi , the integration over (j)i results in merely 
some numerical factor that can be taken care of by a 
proper normalization of the partition function. For the 





/ Of 


9f \ 






at 




df 






di 1 



dl df 
dz dz 



with the derivatives with respect to the Grassmann pa- 
rameters ^ and ^ being understood to be the right ones. 
Evaluating the superdeterminant 



sdetl 



d{ai,aiJ,f) I 



and substituting of 110() into (jSl we are led to the su(2|l) 
path-integral representation of Zt-j given by l@J. Note 
that the U(l) gauge field pi drops out from representation 
An attempt at decoupling the physical electron as 
a U(l) gauge invariant "dressed" holon and spinon has 
been made in Ref. [l9l 

Geometrically, the set (z,^) appears as local (inhomo- 
geneous) coordinates of a point on the supersphere de- 
fined by equation Representation Q-JHI rigorously 
incorporates the local NDO constraint at the apparent 
expense of a more complicated compact phase space for 
the projected electron operators. 



IV. SYMMETRY 

At the supersymmetric point, J = 2t, the t-J model 
Hamiltonian is known to exhibit a global SU(2|1) sym- 
metry. Away from that point this symmetry reduces to 
SU(2)x U(1)C SU(2|1). This symmetry group acts on a 
point (z(t),^(t)) e CP^I^ in a way that. 



Z{t) ^ Zg{t) = 



uz{t) 



—vz[t) + u 



5 6 SU(2) xU(l), 



-vz{t) 



(11) 
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where the group parameters are to be taken to be site 
independent: 



u V 
— V u 



e SU(2), e*^ e U(l). 



(12) 



It can easily be checked that both the SU(2|1) mea- 
sure and the effective action ^ are invariant under 
the group transformations l|ll() . so that the representa- 
tion of the partition function I^J remains intact. No- 
tice that (|11|) appears as a covariant reparametrization 
of CP^''^. However, one can in principle employ any 
other reparametrization, not necessarily of the form of 
the SU(2|1) action on CP'^'^. We are interested in the 
one that decouples the SU(2|1) measure factor into the 
SU(2) spin and the U(l) spinless fermion measures, 

respectively. 

Such a reparametrization can be taken to be 



(13) 



Up to an inessential factor which redefines a chemical 
potential, we get 

DfJ.su{2\l) ^Msn(2)(^,2) X i:'^„(i)(C,C), (14) 
and the effective action becomes 
1— r^- 



St- 



■E 



— I 

1 + ZiZi 



1 



0. 



E / iU^-U)dt- / Hf_jit)dt, (15) 



/3 



with 



J 



5^(i(z,|z,)|2-i)(i-a,0 (i-e^).(i6) 



Here {zi\zj) stands for an inner product of the su(2) co- 
herent states. 



{z.,\Zj) 



1 -t- ZiZj 



From eqs. (|ll|l one can infer the transformation prop- 
erties of the new CP^'^ coordinates (fT^ under a global 
SU(2)xU(l) action: 



, , , , uz(t) + V 

z{t)^zg{t) = _\; 

—vz{t) + u 



(17) 



where 



Note also that 



In 



^ \z)g = e~^'^^\zg). It can be straight- 
forwardly checked that both the measure and the t-J 
action H15|l remain invariant, under such an action of 
SU(2)xU(l). 

The following remarks are needed at this stage. First, 
in spite of the fact that the function {(pg bears a site- 
dependence through the fields, the transformation 
(|17|l is a global one: the group parameters (u, v) are 
site-independent. Second, although the measure factor 
gets decomposed into the su(2) spin and spinless fermion 
pieces, the underlying phase space is not reduced into a 
direct product of the classical spin and a flat fermionic 
phase spaces. The function that enters the transfor- 
mation law for the fcrmions also depends on the spinon 
coordinates, Zi{t). Besides, the symplectic one- form (ki- 
netic term) in the effective action (I15II is not a simple sum 
of purely fermionic and spin contributions. This means 
physically that, in general, the corresponding field exci- 
tations are not independent of each other. In the other 
words, the spin-charge separation does not merely re- 
duce to a simple (z, ^) representation, and should, in 
fact, be described by nonlocal 'string' excitations to be 
constructed out of the basic (z, ^) fields. 



V. EFFECTIVE ACTION 

The spinon amplitudes Zi{t) and the spinless fermion 
parameters, ^i{t), are in fact related to each other by 
the SU(2) transformation laws (|17|l . From this it follows 
that we can construct classical images2£ for the operators 
that describe doped holes. In this respect we make the 
following ansatz: 



*i = 



*t = 



z^ 



zl 



(18) 



It then follows that ^i^*! +l'xvl/^ = = J^', where 
(5^' stands for a classical image of the hole-number oper- 
ator 8 — 1— He = ^ — ^Ylia^"" ■ Therefore the resulting 
fermionic amplitudes describe the propagation of doped 
holes restricted to the NDO constraint. In view of the 
group transformations SU(2)xU(l) for z(t) and the 
^o- amplitudes transform in a linear spinor representa- 
tion of SU(2) as true fermionic amplitudes. Namely, 



u —V 
V u 



(19) 



In terms of the and z amplitudes we get the corre- 
sponding exact representation of the t-J partition func- 
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tion, 



where 



Skin — / H^_j{t)dt. 

Jo 



,(20) 



(21) 



The SU(2) invariant product of the 5 - functions ensures 
the preservation of the correct number of degrees of free- 
dom. The square roots in the S- function arguments come 
from the evaluation of the Jacobian. In this new dopant 
carrier representation the kinetic term, 



Sk 



1 + ZiZ, 



(22) 



is nicely decoupled into purely fermionic and spinon 
parts. It is clear that the fermionic symplectic one- 
form (the first term in (|22|l ) determines a standard 
fermionic symplectic structure d^'o- A d^o- which in 
turn determines the standard Poisson brackets relations 

— 0. As a result the 
corresponding operators ^'^ , ^o-' describe indeed well- 
defined fermionic excitations - in our case, doped holes. 

As a result, using the new fermion fields, the Hamilto- 
nian that corresponds to HfLj takes the form 



(^.^. - 7) 



A' 
,(23) 



where we have dropped the tilda sign. The components of 
the operator of spinon magnetic moment ~§ are the su(2) 
generators in the 5 = 5 representation . Their classical 
images are the components of 5"^' = (2|5'|z) with {S'^)ci = 
|. The hole spin operator A? i '•^t^aa'^a, A? = 
|(5 ^2 - ?j = '^^•^T^ + transforms under ((T7|l . 

as a SU(2) vector while the total hamiltonian (|23|l is a 
SU(2) scalar. It can also be checked that the electron 
spin moment is a linear combination of the above two 
operators: If we integrate out the fields 

in Eq. with the help of the 5— functions, 

we will return to our initial representation as given in 
Eq. 113). 

Different parts of the Hamiltonian H23I) can be associ- 
ated with the different phases of the t-J model. For a 
half-filled band, Eq. (|23|) reduces to the Heisenberg an- 
tiferromagnet (AF), HfJ'j = JJ2^iC^^^j - j)- Away 



from half-filling in the lightly doped regime, where S is 
small enough, so that one can ignore a direct hole-hole 
interaction, the Hamiltonian 



ijij ij 



(24) 



is able to describe the AF - pseudogap (PG) transition 
on a square lattice. Accordingly, at higher doping, the 
Hamiltonian 



ttPG-SC 



(25) 



is appropriate for the pseudogap - superconductor bound- 
ary region of the t-J phase diagram. 

The PG phase itself can be described within our ap- 
proach by a simple SU(2) invariant spinon-fermion inter- 
action. 



(26) 



To see this, one can recast the Hamiltonian (|26ll into 
the form of the phenomenological boson - fermion model 
(BFM) which is known to capture the main observable 
characteristic of the PG phenomenon, namely the reduc- 
tion of the fermionic density of states at the Fermi levelSi. 
Using the Holstein-Primakoff (HP) representation of the 
spin operators on the bipartite lattice L — A^B, 

S^ = 1/2 -bib,, S+^b,, Sr=bl teB, 

where [6,5^] — 1, and performing the following unitary 
transformation of the fermionic operators 

^ *T*, i e A, ^ i e B, 

one is led to the BFM - type Hamiltonian 

^ H^'^^' = t E *i + - J E 

ij i 

- 2jE&Ife. + ^E(*Tj*iA^ + M 

i ij 

+ AE(2fol6. + *l*..-2), (27) 
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\ I r 

t>0 t<0 




FIG. 1: Tc{S) for negative and positive t. The curves from 
the bottom to the top correspond to J/\t\ =0.6, 0.8, 1, 1.2. 



with the impHed summation over a. The Lagrange mul- 
tipher A has been introduced to enforce the constraint 
which assures the vanishing of the total electron spin pro- 
jection, < 2Q^ >= 0. Notice that due to the global U(l) 
invariance the conditions < >— are automatically 
satisfied. While both the conventional BFM and Hf_pj. 
Hamiltonians possess the same global symmetry, namely 
SU(2)x U(l), with the SU{2) group describing of the ro- 
tation of the spinor fields, the origin of the U(l) symme- 
try is different in the two cases. In the standard BFM, the 
global U(l) symmetry corresponds to the conservation of 
the total charge of bosons and fermions. Here the U(l) 
symmetry group appears just as a subgroup of the explic- 
itly broken (by the HP representation) total spin rotation 
group, generated by the operators Qs. Therefore, despite 
the formal similarity between those two effective models, 
the physical contents of the standard BFM and the rep- 
resentation H27|l are indeed different from each other . 

Note, finally, that our discussed hierarchy of effective 
t-J Hamiltonians is basically qualitative in the sense that 
the constraint imposed by the (5-functions might change 
their detailed forms. While this constraint does not seem 
to be of the crucial importance at very low density of 
dopant carriers, it definitely becomes more important as 
S increases, and this may in turn substantially affect the 
final form of the effective interactions. However, since 
the global SU(2) invariance of the t-J Hamiltonian l|23|) 
is not affected by the constraint, it is plausible to sug- 
gest that those changes will at most, at the mean-field 
level, result merely in a renormaization of the Hamilto- 
nian parameters similar to what happens to the t and 
J parameters in the mean-field Plain Vanilla theorjiSS. 
We intend to address these problems in more detail in a 
separate publication. 



VI. SC PHASE 

The RVB mean field treatment of the SC phase of the 
Hamiltonian H23|) is now based on the following assump- 
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FIG. 2: Tc as a function of doping for t < 0. The curves from 
the bottom to the top correspond to J/\t\ =0.35, 0.4, 0.45, 
0.5. For comparison, the insert shows experimental data for 
Naa:Co02 ■ J/H2O taken from Refi-. 

tions: 

i) the global SU(2)xU(l) symmetry is spontaneously 
broken by a local order parameter down to SU(2). The 
SU(2) symmetry is the exact symmetry of the SC phase; 
ii) the dynamics of the SC phase is governed by the BCS- 
type dynamics of the valence bond hole SU(2) singlet 
pairs and is determined by the linearized hole-hole in- 
teraction JX^ii 

as well as by the hop- 
ping term. The hole spin singlets interact with the quasi- 
classical spinon background JJ2ij {^i^j ~ k) 

induced moment-moment interaction J'Yliij ~§ i^j- This 
can be treated within the MF approximation as well. 
However, for the cobaltates we ignore, at first, quan- 
tum fluctuation effects of the spinon field in compari- 
son with the ones originated by the 5*— fieldi^^ In con- 
trast, for the cuprates the {§ iMj) correlation functions 
seem to be of crucial importance and should therefore 
be treated beyond a MF approximation. This is con- 
firmed by the observation of antiferromagnetic ordering 
associated with the superconducting vortex cores^ iii) 
the constrained RVB hole-singlet annihilation operator 
takes the form B^j = ^'it^'jj. - *a^jT' = 0. In the 
SC phase the U{1) global symmetry ^ j„ — + e^^'^ja is 
spontaneously broken by the local SU(2) invariant order 
parameter Ay = (-By ) . 

The t-J Hamiltonian function in Ea. H21|) now reduces 

to 

Hf^j = t ^ + h.c. + - N^^, 

i'ja 

+ Y E (*^T*ji - *a*jT) + h.c. (28) 

where the chemical potential /i has been introduced to 
control the number of electrons, ^ N — ^fjcr. 
This Hamiltonian continues to be invariant under the 
SU{2) action induced by ((T^ . 
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Despite its similar appearance with the standard MF 
RVB result^, ea. (|28|l has a different content. The most 
important difference is related to the presence of con- 
straints imposed by the ^-functions in Eq. H2UI) . How- 
ever, even if these functions are neglected, the Hamilto- 
nian H28(l deals directly with the dopant-particle opera- 
tors ^'io-'s rather than with the electron operators Cia's. 
In terms of ^'^o-'s, N^. has a different representation, and a 
different equation for the chemical potential follows from 
that. As a consequence, within our approach a nonzero A 
directly implies superconductivity in contrast to the orig- 
inal BZA approach, where at half-filling A^^'^ is non- 
zero, but the state is insulating. 

Due to the fact that we are directly dealing with the 
dopant particles vanishes for (5 = 0. Moreover, the 
average of the kinetic term in Hf£^j ^ t5. The increase 
of 5 reflects itself in the gain of kinetic energy which will 
eventually be of the order of J. When this occurs the sin- 
glet pairs tend to break up and the SC phase disappears. 
Strictly speaking Tc is non-zero for any non-vanishing 
5 (see Fig. 1). This is an artifact of the applied MF 
approximation since it neglects a possible onset of mag- 
netic orderingii In the standard RVB decoupling scheme 
the J-term vanishes above the RVB transition temper- 
ature. Therefore, the resulting phase diagram for the 
square lattice usually does not include the AF phase for 
the half-filled case. In the present approach the J-term is 
expressed as a sum of spinon, dopant and spinon-dopant 
terms (see Eq. (I23II '). The introduced above RVB state 
naturally disappears at half-filling and the Hamiltonian 
(|23|l is reduced to the Heisenberg one. Therefore, for 
(5 = 0, the system becomes insulating and the ground 
state energy can only be lowered by the onset of the mag- 
netic phase. In particular, for the square lattice one ex- 
pects the onset of the long range order antiferromagnetic 
order. In order to determine the actual boundary of the 
magnetic phase, one should consider it together with the 
RVB phase, since these phases compete with each other. 

A few remarks are in order. Since the RVB singlets are 
doped hole-hole pairs the present MF favors larger hole 
doping in contrast to the BZA scheme. Thus the present 
approach must be more reliable for the t-J interaction 
on a triangular lattice. 

It can be shown (see Appendix) that the equations for 
the order parameter and the chemical potential that fol- 
low from the partition function representation H20|l with 
the MF BCS Hamiltonian (|28|l are invariant under the 
change t —t, S ^ 1 — S, /i ^ — /i. Thus the NDO 
constraint imposes within the MF BCS approximation a 
symmetry restriction on a possible structure of the phase 
diagram. Namely, the phase diagrams TdS) at t > and 
t < must be located symmetrically with respect to the 
point (5 = i. Any renormalization of the type t — > St, 
frequently used in order to implement the constraint of 
no double occupancy in the MF BCS scheme, evidently 
spoils this symmetry. 

Finally, Eq. (|28|l corresponds to hole doping. How- 
ever, as already mentioned earlier on, the Co02's are 



more likely electron doped compounds. In order to deal 
with this case, within the representation H28|) . we make a 
canonical transformation "^cr '^-a keep the NDO 
constraint as before. In its new form the operator cre- 
ates a dopant electron. The phase diagram Tc x S which 
follows from our new "dual" RVB scheme is shown in 
Fig. 1, for hole doping. If we replace 6 ^ x,t ^ —t we 
reproduce the main figure for the electron doping case. 
Our results for this case are shown in Fig. 2. In the 
insert of this figure we reproduce the experimental data 
from Schaak et al^ for comparison. 

This phase diagram is evaluated directly from Eq. H28() 
considering a triangular lattice of C0O2S. The d + id 
symmetry of the MF OP predicted earlier in Refs. 00 
and I25I is employed throughout the calculations. Other 
symmetries can be tested if necessary using the same 
scheme. The representation (|2U|) with the Hamiltonian 
function given by Eq. H28|) incorporates the NDO con- 
straint rigorously and tells us that at most one spinful 
fermion can live on a given lattice site. Technically, the 
problem reduces to a computation of the fermionic de- 
terminant in the presence of the constraints imposed by 
the ^-functions. 

The fermionic determinant arises upon integrating a 
bilinear form in the exponential over the complex spinors 

=(*r r ). Heren7„ = |(2n+l) 

stands for the Matsubara fermionic frequency and vector 
k G BZ. Had there been no (5-functions in (|20|l the am- 
plitudes S-r!- and E-r^ would have been com- 
k ,CT„>o k ,CT„<o 

pletely independent and contributed equally to the parti- 
tion function. In the presence of the (5-functions, however, 
those amplitudes are no longer independent. 

The (5-functions result in some interference between 
these amplitudes reducing the total contribution to the 
partition function. In order to estimate this reduction at 
the mean-field level we use the following trick. We multi- 
ply the piece of the free energy that comes from the evalu- 
ation of the determinant at the absence of the constraints 
by a coefficient k < 1. Then, requiring that resulting 
equations for the order parameter and the chemical po- 
tential be invariant under the change t — > — t, (5 — s- 1 — (5, 
and /i — !■ —/X gives k = 1/2. 

Although this approximation cannot be justified rigor- 
ously, it goes beyond the one based on the renormaliza- 
tion of the hopping term in the form, t St, which is fre- 
quently used to partly take into account the restriction of 
no double occupancy. In particular, our approximation 
does not spoil the already mentioned symmetry of the 
MF phase diagram under the changes t — > —t, S ^ 1 — S 
dictated by the NDO constraint (see Fig. 1 and the Ap- 
pendix). However, a more detailed analysis must take 
into account a rigorous treatment of the delta-function 
contribution. 

Our results for Naj;Co02 • yH20 are very suggestive 
since the experimentally observed dome structure of the 
phase diagram is reproduced by theoretical calculations 
within a RVB framework. The obtained widths for the 
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dome are also of the same magnitude as given by ex- 
perimentSfSi^ although our doping values are somewhat 
shifted toward the origin. However, recent experimental 
results^^ indicate that the actual hole concentration in 
the cobalt planes may differ from that estimated solely 
on the basis of the No content and the optimal_doping 
can be shifted from the value reported in Ref 2. The 
precise value of J for this compound is still unknown. 
However for t = — O.leV and J/\t\ ranging from 0.35 to 
0.5 as depicted in Figure 2, max varies roughly from 
IK to 4if . Our mean field results are, therefore, in good 
agreement with the existing experimental data. 

The obtained phase diagram is asymmetric with re- 
spect to the change f — > — i (electron and hole doping). 
The t ^ —t asymmetry has also been obtained in Ref. 6 
within a MF slave-boson Hamiltonian. In both MF ap- 
proaches this asymmetry is an obvious consequence of the 
free-particle dispersion relation on the triangular lattice. 
In our case, this asymmetry concerns only the different 
values of the optimal doping in electron- and hole-doped 
systems. In Ref. it is associated predominantly with 
the different width of the SC region in the Tc x S phase 
diagram for different doping regimes. Additionally, we 
have obtained much larger value of the optimal doping 
than that reported in Ref. Note, that the maximal 
value of Tf, obtained in Ref. for the case of electron 
doping is close to that obtained for the hole doping. 

According to our knowledge, there is however no exper- 
imental verification concerning the explicit form of this 
asymmetry in contrast to the electron-hole asymmetry 
observed in cuprate superconductors. The t-J Hamil- 
tonians J^l on a square lattice with double/zero occu- 
pancy for hole/electron doping are unitary equivalent. 
Accordingly, on a square lattice there is no asymmetry 
with respect to the change t —t. In fact, the doping 
asymmetry in the high-Tc cuprates has quite a different 
origin. A possible resolution of this puzzle has rece ntly 
been provided within a two-species t-J model in Ref. [23. 

The doping dependence of the superfluid stiffness Dg 
is an important ingredient of the standard RVB theory 
as discussed in Ref. In particular, small values of Dg 
for 6—^0 determine the superconducting transition tem- 
perature. This quantity can be obtained with the help of 
the linear response theory from the relation between the 
current and the transverse gauge field. The response ker- 
nel consists of paramagnetic and diamagnetic parts. In 
the superconducting state (or more precisely for A 7^ 0) 
the paramagnetic contribution vanishes for T 0^22, It 
has been shown that the diamagnetic part imposes also 
the upper bound on Dg-^^ In the case of the hypercu- 
bic lattice with the nearest-neighbor hopping, Dg/n is 
bounded by the absolute value of the kinetic energy on a 
bond, whereas for a more general dispersion relation e(k) 
the kinetic energy should be replaced bji^S 
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FIG. 3: Doping dependence of K{T) and A(T) obtained for 
the triangular lattice with J/\t\ —0.5 and T — + 0. 



where m^^{k) = d^eik) / dkadkjj. In order to estimate 
the magnitude of the superfluid stiffness in our approach, 
we have calculated T = limit of this quantity. Fig 3. 
shows the doping dependences of K{0) and A(0). In 
the standard RVB approach the NDO constraint t ^ 5t 
renormalizes the superfluid stiffness since it modifies the 
kinetic energy. Consequently, close to half-filling Dg 
vanishes despite the finite value of AiSS In this regime 
Tc Dg S. In our approach the superfiuid stiffness 
vanishes for the half-filling as well. Since the kinetic en- 
ergy term in Eq. H23|l contains only the dopant particles 
the vanishing of the superfluid stiffness for (5 —> is an in- 
trinsic feature of our approach and occurs independently 
of the applied approximations. 

We end this section discussing briefly how one can 
control the BCS mean-field decoupling H28|l within the 
su(2|l) supersymmetric representation of the Hubbard 
operators. This can be done by means of a large-iV ex- 
pansion based on a generalization of the SU(2) globally 
invariant t~J Hamiltonian ^ in terms of the symplec- 
tic group Sp(2N) of 2Nx 2N unitary matrices (note that 
Sp(2)^ SU(2))iSi2L2i: 

+:^E^-'^'7.'^'^'''^"'''-A^E^r, (29) 

ij i 

where the summation over the Sp(2N) indexes a,r] — 
±1, ±2, zLN is assumed and the Sp(2N) antisymmetric 
tensor ecrcr'= sgn(a)6cr,-a' ■ The local constraint of the t- 
J model can now be taken in the form V Xf"' — N at 
every lattice site. 

The exchange term in (|29(l can be rewritten in terms 
of the Sp(2N) invariant valence - bond operators B^^ := 

Y.aa' eaa'Xfxf"" in thc form 



K{T) 



I 

TV 



n(fc)Tr 

keBZ 
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and may be decoupled by the link field 

aa' 

Ki N — CO this decoupling becomes exact. 

The su(2|l) supersymmetric representation of the Hub- 
bard operators can be extended to the case of the 
su(2iV|l) superalgebra. Since Sp(2N)c SU(2N) one can 
then employ the su(27V|l) coherent states and the corre- 
sponding path integral to treat the t-J Hamiltonian (|29|l . 
In this way we will eventually arrive at an Sp(2N) glob- 
ally invariant generalization of the representation (|23|l . 
with S and M operators being now replaced by the cor- 
responding Sp(2N) generators. More details will be given 
elsewhere. 



VII. CONCLUSION 

To conclude we developed a RVB mean field the- 
ory which takes a direct account of the dopant carri- 
ers. These dopant particles are represented by appropri- 
ate fermion fields which carry both spin and charge and 
transform themselves as true SU(2) spinors. The result- 
ing theory is written in a very convenient form since we 
are able in this way to consider the both doping depen- 
dence of the critical temperature as well as the kinetic 
energy effects which eventually destroy the superconduc- 
tivity at larger dopings. By making a more extensive 
use of Hubbard operators we go beyond the conventional 
slave-boson approximation and take sufficient care of all 
symmetry properties of the Hamiltonian model. Since 
we apply a mean-field decoupling, there is always a gap 
in the energy spectrum of the dopant particles when- 
ever A 7^ 0. Consequently, the superconducting tran- 
sition temperature kT^ and the energy gap JA are of 
the same order. In the slave-boson RVB formulation of 
Ref. 0, for low doping corresponds to the condensation 
of holons. Therefore, in that approach is decoupled 
from the value of A. Such a decoupling may also occur 
within our approach, e.g., when spinous are considered 
beyond the mean-field level. One can see from Eq. H30|) 
that the lowest-order spinon-holon coupling takes on the 
form ^ Yliij CiCj i^j — Zi) + h.c. The resulting density of 
states for holons may be finite also for A ^ Oi^ 

We initially applied this new RVB scheme to describe 
the superconducting properties observed in the cobal- 
tates. We succeeded in getting qualitative good agree- 
ment with experiment. The dome structure of the phase 
diagram Tc x S is well reproduced within a RVB frame- 
work33.. This is achieved without any symmetry violation 
of the t-J model for the whole doping regime. 

While preparing this version of our work we came 
across another RVB formulation in terms of dopant 
carriers^. Those authors use an extended t-J model with 
t, t' and t" hopping parameters. In their scheme how- 
ever those parameters are renormalized by interactions 



and this procedure automatically violates the underlying 
symmetries of the original t-J model. 

As discussed in the Appendix, the MF phase diagram 
Tc{S) for the t-J model on the square lattice without 
frustration results in a max Tc located a.t 6 = 1/2. Al- 
though incorporating the next-nearest-neighbor (NNN) 
interaction in the kinetic term slightly shifts the diagram 
toward the origin, it cannot account for the experimen- 
tally observed curve for the cuprates. It's clear that frus- 
tration is an important ingredient for the success of our 
RVB method. However, apart from that there is yet an- 
other important feature with needs to be taken into con- 
sideration to properly deal with the cuprates case. In 
the cuprates there are strong antiferromagnetic correla- 
tions which manifest themselves even inside the super- 
conducting vortex cores. As a result the < '^'iMj > 
correlations, which seem unimportant for the cobaltates, 
may also play an important role in the cuprates. This 
will produce strong phase fiuctuations which, most likely, 
need to be taken into account beyond mean field approx- 
imation. This work is in progress and will be presented 
elsewhere. 

As a final concluding remark let us say a few words 
about the ideology of the present paper. The basic idea 
is to use Hubbard operators, instead of the standard 
fermion operators accompanied with the nonholonomic 
constraint of no double occupancy. This enables us to 
impose the NDO constraint locally at each lattice site. 
This constraint results in strong electron correlation ef- 
fects which are believed to be essential ingredients for 
dopped Mott insulators. Since the Hubbard operators 
appear as Gutzwiller projected (GP) electron operators 
on the states with no double occupancy, it is in princi- 
ple reasonable to work directly with the GP operators 
and wave functions. In this way Paramekanti, Randeria 
and Trivedi recently studied the Hubbard model making 
use of parameters relevant for the cuprates, in the frame- 
work of the variational Monter Carlo GP d-wave stated. 
They showed that the strong electron correlations im- 
posed by the Gutzwiller projection destroy the off di- 
agonal long range order as S ^ qualitatively tracking 
the observed nonmonotonic Tc{S). Basically the same re- 
sult follows from the Plain Vanilla version of RVB, where 
the Gutzwiller projection is treated within the mean-field 
representation. 

In our approach we also treat the NDO constraint at 
the mean-field level. However, we go a step further since 
we make explicit use of the algebraic relations between 
the Hubbard operators, namely those of the su(2|l) su- 
peralgebra. This adds some extra information which 
is encoded in the superalgebra commutation relations. 
In particular the classical phase space realization (the 
coherent-state representation) of su(2|l) provides us with 
the complex canonical coordinates (z, ^) which eventu- 
ally appear as the basic spinon-fermion fields in the path 
integral effective action (0}. Dopant quasiparticle ampli- 
tudes (|18f) are constructed out of these fields appearing in 
our theory as emergent phenomena. We arrive naturally 
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in this way at the RVB theory for dopant carriers. 



Acknowledgments 

One of us (E.K.) wants to acknowledge the hospitahty 
of the ICCMP's staff and the financial support received 
from CAPES - Brazil. 



VIII. APPENDIX 

In this Appendix we show that within the BCS MF 
approximation H28|l the equations for the order param- 
eter and the chemical potential are invariant under the 
change t —t, S 1 ~ 6, /i —fi, provided the NDO 
constraint is rigorously taken into account. 

First, we integrate out the fields in Eq. 

with the Hamiltonian function given by (|28|l . which re- 
sults in the effective action (|15|) with the classical Hamil- 
tonian function now being, 



TTCl 

^sc 



JA 



-t^(6,Cj(^»kj) + h.c.) ~ 



2 5 VOTRFKTT^ 



+ h.c] . (30) 



Here Zi{t) and ^i{t) are dynamical fields. This represen- 
tation rigorously incorporates the constraint of no dou- 
ble occupancy. Because of the rather complicated form 
of the action, we are in general unable to write out ex- 
plicitly a quantum counterpart of Hamiltonian l|3U|) as a 
function of the su(2) spin generators and spinless U(l) 
fermionic operators. However, in the SC phase we get 
Zi{t) — Zi, which means that quantum fluctuations of the 
background spinon fields are ignored. In that case only 
the fermionic kinetic term is left in the action (|15|l . and 
the quantum Hamiltonian can be easily identified, 

Hsc = -t + h.c.) - /i ( iV - ^ J 

ij \ i / 



J A Iff Zj — Zi 



+ h.c 



(31) 



The fiS stand for the on-site spinless fermionic operators, 
with {/i, /J} = 6ij that correspond to the classical Grass- 
mann amplitudes, —: ^, which give = 0. The 

dynamical spinon field Zi(t) looses its time-dependence 
and turns itself therefore into a sort of external classical 
c- valued spinon field. 

Next, we evaluate the on-site free energy function. 



F/N 



N 



(32) 



where the symbol Tr is used to indicate the summation 
over the fermionic degrees of freedom as well as the com- 
plex c- valued spinon fields: 



Tr {■■■):= j Dfisu{2){z,z)trfjf {■■■) (33) 

The z-integral in (|33|l appears as an ordinary multiple 
integral. In this way the order parameter and chemi- 
cal potential are determined by the conditions dF/ dA = 
0, dF/d^i = 5 — \ which explicitly give 



Zj Zi 



+ h.c) = 0, (34) 



and 



(35) 



respectively Here ((• • •)) := Tr {■ ■ .)e-'^"sc jXr e-f^""^^ . 
It can be checked straightforwardly that eqs. H34|l and 
(|35|l are invariant under the change t — > — ^ — > 
— /i, (5 ^ 1 — (5. To see this one should simultaneously 
make the canonical transformation, fi fj, and change 
the integration variables, Zi — > —Zi. Accordingly, the 
phase diagrams Tc{5) at < > and t < are located 
symmetrically with respect to the point S = 1/2. 

Explicitly, the equations for the order parameter and 
chemical potential read 



1 

N 



E 

keBZ 



tanh(^) 



1 ^ tanh(^)^ ^ ^ 



z 
J' 



keBZ 



Ek 



where 



Ei^itk-i^r + J'^'\Pk 



(36) 



(37) 



(38) 



and — -~'2tj^, — ^-cos/cn. In the case of the 2D 
square lattice 7^ = cos + cos ky , whereas 7^ = cos kx + 
2 cos{kx/2) cos{ky\^/2) for the 2D triangular lattice. For 
the dx2-y2 pairing on the square lattice the phase factor 
reads Pj: — coskx — cosky. For the triangular lattice we 
assume a. di + id2 symmetry of the order parameter^i^ 
Then, 



k n,. 

COS — cos ^ 
2 



fc.„V3 



kx . ky\'^ 



iv3 sin — sin 

2 2 



/Sr = COS kx 

(39) 

The equations (|36|l and (|37|) are clearly seen to be invari- 
ant under the change t —t, — > — /i, (5 — > 1 — (5, which 
results in the phase diagram depicted on Fig.l. 
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Note that the t-J Hamihonian on a square lattice with 
the nearest-neighbor (NN) interaction is invariant under 
the change, t ^ —t. This is because this change amounts 
to a certain unitary transformation of the lattice electron 
operators. It then follows that the above two phase di- 
agrams merge in this case into one, located at (5 = 1/2. 
Incorporating frustration (e.g., by taking into account 
the NNN interaction in the t-dependent term) destroys 
this symmetry and results in splitting of this diagram 
again into two located symmetrically with respect to the 
point S = 1/2. However for the generic values of the t-J 
parameters that splitting is rather small and cannot ac- 
count for an experimentally observed phase diagram for 
the cuprates. 

If we ignored completely the NDO constraint taking 
into account the modes S-^ and S-r> on equal 

grounds, we would get (on a square lattice) a diagram 
with max Tc located at S = I. This is markedly different 
from the NDO constraint-free BZA result, where max 
occurs at (5 = 0, which bears out that our theory is in a 
sense dual to the original BZA approach. 

The conventional BZA MF theory formulated in terms 
of the lattice electron spin singlets with the renormaliza- 
tion t St being implemented to partly incorporate the 
NDO constraint, however fails to maintain the symmetry 
of the phase diagram dictated by this constraint, and re- 
sults in the same observation: max Tc takes place again 



at (5 = 0, as in the constraint-free BZA theory. To see 
this consider the BZA MF Hamiltonian,^ 

ijG i<7 

+ — ^ (c,TC,i - cacjt -I- h.c.) + ^.(40) 

One obtains the following system of equations to deter- 
mine the order parameter and chemical potential: 

Et ^fe 2J' ^ ' 

k " 

k 

where 

i?| = (i,,-M)' + j2A^|, (43) 

and tj: = —2tSj^, = ^-cosfcn. The ensuing phase 
diagram T^^^{6) is invariant under the change S — > —S 
so that max Tc always occurs at 6 — 0. 
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